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ABSTRACT 

We use X-ray morphological measurements and the scatter of clusters about observed and simulated 
scaling relations to examine the impact of merging and core-related phenomena on the structure of 
galaxy clusters. We use a range of X-ray and near-infrared scaling relations; all observed scaling 
relations constructed from emission-weighted mean temperature and intracluster medium mass, X- 
ray luminosity, isophotal size, and near-IR luminosity show a separation between clusters identified 
as cool core (CC) and those identified as non-cool core (NCC). We attribute this partially to a simple 
temperature bias in CC clusters, and partially to other cool core-related structural changes. Scaling 
relations constructed from observables that are largely independent of core structure show smaller 
separation between CC and NCC populations. We attempt to minimize cool core- related separation 
in scaling relations via two methods: by applying a uniform scale factor to CC cluster temperatures 
and determining the scale factor for each relation that minimizes the separation between CC and NCC 
populations, and by introducing cluster central surface brightness as a third parameter in observable- 
temperature scaling relations. We find an average temperature bias factor of 1.07 ± 0.02 between the 
CC and NCC populations; the three parameter approach reduces scatter in scaling relations more 
than a simple CC temperature scaling. 

We examine the scatter about the best-fit observable-temperature-brightness scaling relations, and 
compare the intrinsic scatter within subsamples split by CC/NCC and four different morphological 
merger indicators. CC clusters and clusters with less substructure generally exhibit higher scatter 
about scaling relations. The larger structural variations in CC clusters are present well outside the 
core, suggesting that a process more global than core radiative instability is at work. Simulations 
without cooling mechanisms also show no correlation between substructure and larger scatter about 
scaling relations, indicating that any merger-related scatter increases are subtle. Taken together, 
the observational and simulation results indicate that cool core related phenomena — not merging 
processes — are the primary contributor to scatter in scaling relations. Our analysis does not appear 
to support the scenario in which clusters evolve cool cores over time unless they experience major 
mergers. 

Subject headings: galaxies: clusters: general — X-rays: galaxies: clusters — intergalactic medium 



1. INTRODUCTION 

Galaxy clusters provide a setting for exploring the com- 
position of the universe on large scales, and for studying 
the growth of structure. X-ray and Sunyaev-Zel'dovich 
effect surveys of large numbers of clusters will soon be 
used to study the nature of the dark energy. It is thus 
vitally important to understand cluster structure and its 
connections to observable, bulk properties of clusters. 

Clusters exhibit strikingly regular power law scaling 
relations between such properties as emission-weighted 
mean temperature, intracluster medium (ICM) mass, 



IR light (e.g., 


Mohr & Evrardlll997l IMohr et al.lll999l: 


IVikhlinin et alJ 


2002HLin et al.ll2003l). For relaxed, iso- 



lated systems this is expected, although even for isolated 
systems there should be some scatter at a given average 
temperature due to variations in formation epoch, gas 
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fraction, and galaxy and star formation history. 

However, clusters are in fact young, dynamic systems. 
There is evidence, much of it based on quantitative sub- 
structure measures, for merging in a significant frac- 
tion ( >50%) of nearby clus t ers (e . g.. iGeller fc Beers! 
T98l IDressler ShectmTl ITgM IMohr et a,lJ 1199,'^ 



Buote fc Tsail 1199(a) . The advent of high-resolution X- 



ray instruments such as Chandra has made observation of 
the complex hydrodynamical structure of merging clus- 
ters possible, revealing and permitting detailed measure- 
ments of ICM properties aro und features such as "cold 
fronts" (Vikh hnin eraI]l200H) . 

This dual nature of clusters — frequent merging on the 
one hand, and regular scaling relations on the other — is 
puzzling. The extreme energetics of merger events, in 
which ~ 10^^ erg of kinetic energy can be thermalized 
and cluster structure greatly disrupted, would suggest 
that cluster properties should not be correlated in such a 
simple way. One might expect to find a statistical corre- 
lation between the deviation of a particular cluster from 
a scaling relation and the substructure — an indication of 
merger activity — in that cluster. 

Tight scaling relations are observed in simulations 
of clusters that ev olve in a cosmological context (e.g., 
lEvrard et alJll996l: iBrvan k Normanlll998l: iBialek et all 
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1200 ll) . However, simulations of mergers of isolated, re- 
laxed clusters suggest that merging clusters should ex- 
hibit massive boosts in cluster bulk parameter s such 
as temperature and X-ray luminosity ( Ricker & S arazinI 
120011 iRandall et al. "2002n. While the simulations indi- 
cate that these boosts are correlated, it seems likely that 
a cluster involved in a major merger would stand out 
from scaling relations constructed from several observ- 
ables simultaneously. 

If the cool cores found in a large fraction of the clus- 
ter population are an expected outcome of clu ster relax- 
ation in the absence of merger activity (e.g., lOta et all 
l2005j) ■ then we may expect "relaxed" clusters to dis- 
play less structural variation than clusters recently in- 
volved in mergers. If, on the other hand, cool core 
structure is linked to more fundamental properties of 
individual clusters, suc h as early-time entropy injection 
ijMcCarthv et al.lf2004|) . then merger-related effects may 
be overshadowed. There is presently no consensus as to 
whether merger-related effects dominate the scatter in 
scaling relations (e.g.. iSm ith et al. 2005)jOr are a mi- 
nority contribution (e.g.. lBalogh et al.ll200^ . 

In this paper we study multiple scaling relations. One 
may envision surfaces of clusters in a hyperplane con- 
structed from several observables, and imagine studying 
deviations from these relations in multiple dimensions. 
We begin by examining only two parameters at a time, 
largely due to the greater ease of visualizing and under- 
standing the deviations. Later we combine these two- 
observable analyses to examine cluster departures from 
the population in the much higher-dimensionality space. 
Our hope is that by examining nine different X-ray and 
near-infrared observables, each of which represents a dif- 
ferent integral over the structure of the ICM and galaxies, 
we will be able to discern even modest structural devia- 
tions and probe their relationship to merging and other 
cluster phenomena. 

If, as expected, cluster mergers perturb crude 
observables — for example, producing erroneously high 
estimates of cluster mass — then this must be taken into 
acco unt in surveys using clusters to study c osmology 
(e.g. Haiman et al."2001; Hu '20031 iMaiumdar fc Mohii 
12003. .2004 .Lima fc Hu 2005) . It is thus of great im- 
portance to understand the true effects of merging on 
cluster bulk properties. Positions on scaling relations 
have been examined for a fe w ind ividu al merging clus- 
ters, e. g., CI J0152.7-1357 p ^ ghan c TaP 120031) and 
A2319 l|0'Hara et ani2004f) . with no significant devia- 
tions found. Clearly, however, analysis of a larger sam- 
ple of clusters is needed to make meaningful statements 
about the relationship between merger signatures and 
bulk properties. Such an undertaking is the subject of 
this paper. 

Scaling relations have already been used to study clus- 
ter structure; for example, by examining the the ef- 
fects of structure on the slope of the relations (e.g., 
Evrard et al. 1996; Cav aliere et a l. 1997: ' Mohr fc Evrardl 
1997: ^Bialek .et aL ■200lE The trend in recent X-ray 
studies of clusters, however, has been toward detailed 
studies of clusters using high-resolution instruments such 
as C handra, and correspond i ngly detailed simulations 
(e.g., i Ricker fc SarazinI 120011: iNagai fc K ravtsov' 120031: 
lOnuora e_Lal,. 2_0_0_3i) . These studies have uncovered many 
surprising aspects of cluster structure. There is still a 



need, though, for a clearer picture of the effects of merg- 
ing on the entire cluster population. An analysis of pop- 
ulations of simulated clusters suggests that it will not be 
possible to isolate "undisturbed" or "relaxed" clusters in 
large samples over a range of redshifts, and so a better 
understanding of merger effects on cluster bulk proper- 
ties and on the general population is required. 

We present here a joint analysis of a flux-limited sam- 
ple of 45 nearby clusters observed with the ROSATPSPG 
and the Two Micron AU Sky Survey (2MASS), and an 
ensemble of 45 simulated clusters evolved in a cosmolog- 
ical context. We begin with the scaling relations from 
the observational sample in § 13 In § 13 we discuss the 
effects of cool cores in our sample. We discuss a method 
of correcting for cool core effects so that we can exam- 
ine only structural, merger-related scatter. We discuss 
the use of peak surface brightness as an indication of 
cool core strength in § 0] § [S] reports the study of the 
relationship between displacement from scaling relations 
and substructure in both observed and simulated clus- 
ters. Finally, we list our conclusions in §|S1 

Throughout the paper we assume a ACDM cosmology 
with r^M = 0.3 and 17a = 0.7, and use a Hubble param- 
eter of Hq = 7O/170 km s~^ Mpc~^. All uncertainties are 
68% confidence, or 1 a. 

2. OBSERVED SCALING RELATIONS 

We study an ensemble of 45 members of the [Edge et aj 
1^990) flux-limited sample, observed with ROSATPSPC. 
We use the same reduced imaging data as 'Kdohr et a]} 
fl999) (hereafter MME). These data have a pixel scale 
of 14". 947, and an energy range of 0.5-2.0 keV. The reso- 
lution of PSPC is, of course, poorer than the current gen- 
eration of X-ray instruments, but none of the observables 
we are measuring require higher resolution. For details 
of the reduction, see MME. For 34 of these c lusters we 
also u se measurements of K-hand light from iLin et all 
l)2004l) : the reduction, measurements, and uncertainties 
are discussed therein. 

We use previously published emission-weighted mean 
temperatures measured with Ginga and ASCA, except 
for A2244, where we use a temperature measured with 
the Einstein MPC, and A3532, for which no appropriate 
temperature with uncertainties can be found in the litera- 
ture, and for which we use a temperature derived from a 
lumin osity-temperature relation (Rciprich fc B6hringe3 
I200I . We use PSF-corrected, cluster X-ray peak sur- 
face brightness values from MME, but all other X-ray 
observable quantities are measured as part of this anal- 
ysis. We divide the sample into cool core (CC) and 
non-cool core (N CC) clusters based on pubhshed cen- 
tral cooling times l|P eres et aLl |199 8'l. We adopt the clas- 
sification of iMohr^TEvrardl l)1997|l in which CC clus- 
ters are those with central cooling times at least 3a be- 
low 7.1 Gyr (for Hq = 70 km s'^ Mpc'^; 10 Gyr for 
Hq — 50 km s^^ Mpc^^). By this measure our sam- 
ple of 45 clusters contains 30 CC and 15 NCC clusters. 
Basic information about the sample, plus the measured 
observables, is given in Tabled 

In this section we first examine the X-ray luminosity- 
temperature relation in detail as an example of the 
scaling relations we are using. We point out features 
that will be common to all relations and demonstrate 
that the source of scatter is true structural variation 
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TABLE 1 
PSPC Sample Information 
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among clusters, not measurement uncertainties. We then 
present the remaining X-ray scaling relations and quan- 
tify their scatter as well. Finally, we present the near- 
IR luminosity-temperature relation and discuss the addi- 
tional information available from the galaxy properties. 

2.1. X-ray Luminosity-Temperature Relation 

We begin by examining the scaling relation constructed 
from Lx{<r5oo) (i.e., the X-ray luminosity projected 
within a distance corresponding to the virial radius rsoo) 
and the emission-weighted mean temperature Tx (all 
temperature values given in this paper are emission- 
weighted X-ray temperatures). The luminosity is mea- 
sured in the 0.5-2 keV band. The virial radius taIs 



the radius within which the mean density is A times 
the critical density of the universe; these radii are ob- 
tained from observed Mvir-Tx virial scaling relations. 
For rsoo we use a mass-temperature relation obtained 
using clusters with mas ses greater than 3.6 x 10^'^ Mq 
l)Finoguenov et al]l200H) . which gives 

rgoo = 0.447 h^^' ^^x''^ Mpc . (1) 

Note that in our analysis we use the emission- 
weighted mean temperature rather than cool core- 
corrected emissio n-weighted mean temperatures. The 
iFinoguenov et all pOOl) relation is derived using tem- 
perature profiles measured with ASCA, but their 
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emission-weighted mean temperatures have been cor- 
rected for cool core effects. 

One CC cluster in our sample, A3526, has a value of 
r5oo large enough that it extends beyond the edges of 
the PSPC image. Three other CC clusters (Ophiuchus, 
A426, and A1060) have values of r^oo that become too 
large for the images when the temperature is scaled to 
account for cool core effects in the analysis presented in 
§ 13 We therefore exclude these four CC clusters from 
our analysis involving ix(<^50o) (and, below, from the 
analysis of the core-subtracted luminosity within rsoo). 

We consider X-ray luminosity uncertainties from three 
sources. First, we measure luminosity projected within a 
virial radius determined by the temperature, so there is 
an uncertainty in the radius due to the uncertainty ATx 
in the measured cluster temperature Tx- We thus mea- 
sure the luminosity projected within radii determined 
from temperatures Tx + ATx and Tx — ATx, and av- 
erage the absolute difference between these luminosities 
and the luminosity measured at the virial radius. There 
is also some uncertainty in the X-ray images; we use er- 
ror images created as described in MME, and sum the 
error within the region of interest just as we do for the 
X-ray image counts. Finally, there is an uncertainty in 
the background, which we account for by raising and 
lowering the background level by 10%, measuring the 
luminosity at each level, and averaging the deviations 
from the value measured with the standard background 
level. The temperature-induced uncertainty. X-ray im- 
age uncertainty, and background uncertainty are added 
in quadrature to obtain the total uncertainty in the X-ray 
luminosity. 
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Fig. 1. — X-ray luminosity projected within rsoo (iop), and 
deviation of clusters from the best fit relation (bottom). Open and 
filled markers correspond to CC and NCC clusters, respectively. 
The uncertainties include both the measurement uncertainty in 
the luminosity and an effective temperature contribution to the 
luminosity uncertainty as described in the text. 

When we examine the deviation of a cluster from a 
scaling relation, we do so in one dimension only, e.g., 
only for luminosity in the luminosity-temperature rela- 
tion. Wc thus include a temperature component to the 
uncertainty in each observable, which we approximate 
using the best fit power law to the observed scaling re- 
lation. That is, for each observable O we determine the 



scaling relation exponent a (i.e., O cx T^), and find a 
fractional uncertainty due to the temperature. 



which is added in quadrature to the measurement uncer- 
tainty. 

Figure^shows the Tx(<?'50o)~'7x scaling relation, in- 
cluding uncertainties. One immediately notices that 
there is a separation of the CC and NCC populations 
on the relation, i.e., the CC clusters tend to lie above 
the best-fit line, and the NCC clusters li e below; this 
relationship has been observed before fe.g.. "Fabian et all 
©94; Markevi tch 1998; McCarthv et al 2004) . One part 
of the separation is likely a simple temperature bias; that 
is, CC clusters have a cool central region that does not 
significantly affect the structure of the cluster outside 
that region, but whose relatively high density leads to 
a significant bias toward lower emission-weighted mean 
temperatures relative to clusters without a cool core. 
The remainder of the separation is due primarily to the 
higher central gas density in CC clusters, which leads to 
greater X-ray luminosity at a given temperature. 

Because we will be examining the scatter about scal- 
ing relations and drawing conclusions about the effects 
of substructure on those relations, it is important to es- 
tablish that the scatter is a real, intrinsic effect, and not 
due to measurement uncertainties. This is made plain 
qualitatively from the bottom portion of Figure ^ which 
shows the deviation from the best-fit relation; the un- 
certainties in Lx are clearly smaller than the intrinsic 
scatter in the relation. To address the issue quantita- 
tively, we measure both the raw (i.e., RMS) scatter and 
the intrinsic scatter, which we express in terms of InO, 
by finding the values of a for which the reduced value 
for a given relation is unity. We find the intrinsic scat- 
ter tJint by adding a uniform value in quadrature to the 
measurement uncertainty for each cluster. 

For the Tx(<?'5oo)-Tx relation, we find CTint^fraw^ 
0.53 (these values are also listed in TableEJ; that is, the 
intrinsic and raw scatter are the same to the precision 
given here, and so the intrinsic scatter is clearly much 
greater than the measurement uncertainties. We may 
thus be certain that the scatter about this relation is 
due to real structural differences between the clusters. 

2.2. Other X-ray Scaling Relations 

In addition to the luminosity projected within rsoo, 
we measure Tx(<?'250o)i the luminosity projected within 
''2500; studying relations at different radii gives us 
information about the effects of structural variations 
on different scales within the cluster. To determine 
^"2500 we use an M2500-TX relation derived from Chan- 
dra observations o f relaxed intermediate-redshift clusters 
l|Allen et al.ll2nOH) . which gives 

r25oo = 0.227 h^^ T^^°^ Mpc . (3) 

We also measure the core-subtracted X-ray luminos- 
ity (Txcs) within rsoo- We exclude the luminosity 
projected within a radius corresponding to 0.20r5oo, 
which minimizes t he effects of core structure (e.g., 
iNeumann fc Arnaud.1999) . This radius corresponds to 
0.13r2O0j where the virial radius r2oo is calculated from 
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TABLE 2 

Raw and Intrinsic Scatter in Scaling Relations 



Observations Simulations'' 

Original With (JO 3 Parameter 

Relations Temp. Scaling (C-Tx-/o) 

Scaling Relation (Jraw CT~t o"raw o\^t ^"t o-int 



MiCM(<r-500)-Tx 


0.19 


0.17 


0.18 


0.17 


0.16 


0.15 


0.20 


Micm(<»-2500)-Tx 


0.23 


0.22 


0.20 


0.19 


0.16 


0.14 




I'x(<?-50o)-Tx 


0.53 


0.53 


0.34 


0.33 


0.26 


0.24 




-t'x(<'"2500)-Tx 


0.69 


0.67 


0.39 


0.39 


0.36 


0.34 




-i'XCs(<''500)-T'x 


0.29 


0.28 


0.25 


0.23 


0.23 


0.21 


0.27 




0.14 


0.14 


0.14 


0.14 


0.13 


0.13 


0.10 


^1.5xlO-13"^X 


0.17 


0.17 


0.17 


0.17 


0.14 


0.14 




i'NIR(<»'500)-T'x 


0.19 


0.19 


0.17 


0.17 


0.19 


0.19 




-^3 X 10- !•* ^-'^ICM (<^500 ) 


0.06 


0.06 












-^3x lO-l^^-'^XCS (<'~500) 


0.06 


0.05 












-^3 X 10- --^'NIR (<»"500 ) 


0.16 


0.16 













Tor simulations, A/icm(<^50o) is actual gas mass, not measured from mock observations, and 
the isophote for Rj is chosen to produce normalization similar to the observed i?3^j^Q-i4-Tx 
relation. 



''soo by using an NFW dark m atter density profile with 
concentration parameter c = 5 ijNavarro et al.lir997j) : the 
relationship is r2oo = l-Slrsoo- 



s 10 




2 4 6 810 
k,T, [keV] 



Fig. 2. — From top to bottom, we show scaling relations for 
A^iCM within rsoo, isophotal size for an isophote of 3 X lO"^** 
erg s~^ cm~^ arcmin"^, the X-ray core subtracted luminosity 
projected within rsoo, and the A'-band luminosity projected within 
''500 1 versus Tx. The open and filled markers correspond to CC 
and NCC clusters, respectively. 



We measure the ICM mass within r^oo f^nd r25oo ■ The 
masses are measured as discussed in MME, using the 
/3-model parameters given in that paper (some clusters 
are fit with a double /3- model), the emission- weighted 
mean temperature, and a measurement of the X-ray flux. 
We use the Micm(<?'50o) uncertainties from MME as 
a starting point, but we adjust the temperature uncer- 
tainty contribution to reflect the newer and more accu- 
rate temperatures available for some clusters. We adopt 
the same fractional uncertainty for Micm(<''250o) for 
-^^icm(<'"50o); the median fractional uncertainty is ~3%. 

We determine the isophotal size Rj of a cluster by mea- 
suring the area Aj enclosed by a particular isophote / 
and finding the effective radius given by Aj = nRj. We 
measure Rj for two isophotes: 3 x 10"^'* erg s~^ cm~^ 
arcmin"^ and 1.5 x 10~^^ erg s~^ cm~^ arcmin"^, in 
the 0.5-2 keV band. The lower isophote lies well outside 
the core region of the clusters, and so, hke the lumi- 
nosities and masses measured within rsoo, reflects clus- 
ter structure in a way largely unaffected by core sub- 
structure. The higher isophote is more reflective of core 
structure. For some clusters in our sample the central 
surface brightness does not rise (or barely rises) above 
the brighter isophote; these eight clusters are excluded 
from all analysis at this isophote. 

When measuring the isophotal size we include a back- 
ground uncertainty of 10%, calculate Rj with the higher 
and lower uncertainty, and take the average of the devi- 
ations from the standard background value to obtain the 
uncertainty in Rj. PSPC also has a ~10% uncertainty 
in its effective area. However, changing the effective area 
tends to simply shift the entire cluster population up or 
down in isophotal size, and does not affect the deviation 
of individual clusters from the relations; hence, we ignore 
this particular uncertainty. In the soft X-ray band the 
conversion from PSPC counts s^^ to physical flux is ap- 
proximately independent of cluster temperature, and so 
temperature uncertainties do not lead to uncertainties in 
the measured isophotal size. 

Each observable derives from a different integral over 
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cluster structure: 



1 



where A(r) is the emissivity of the ICM gas, and the 
other symbols have their usual meanings. Isophotal size 
derives from cluster structure in a somewhat more com- 
plex way; see lMohr et alJ l)2000() for a discussion. Thus 
by studying multiple observables, we are not simply in- 
creasing the size of our parameter space; we are, in fact, 
looking at several different ways of quantifying the struc- 
ture of clusters. By focusing on crude observables like 
these we are able to work in a regime where the mea- 
surement uncertainties are small compared to the intrin- 
sic scatter. 

The scaling relations for AficM(<''50o) ^sxio-^'' ^-^d 
ixcs(<''50o) are shown in Figure |21 in Figure |31 we 
plot relations for Mjcm within r2500j ix(<?'250o)j and 
i?i 5x10-13- For clarity we do not show measurement un- 
certainties; however, as shown for Lx(<?'50o) in Figure^ 
they are significantly smaller than the intrinsic scatter. 
The actual measurements of raw and intrinsic scatter 
are given for all relations in Tabled the intrinsic scatter 
dominates the total scatter in all cases. 
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Fig. 3. — Scaling relations for AfjcM within r2500, X-ray lu- 
minosity projected within r2500, and isophotal size for an isophote 
of 1.5 X 10~ erg s~^ cm~^ arcmin"^, versus Tx. The open and 
filled markers correspond to CC and NCC clusters, respectively. 



All X-ray relations indicate a separation between CC 
and NCC clusters. The non-core-subtracted luminosi- 
ties show the largest effect; this is not surprising, be- 
cause the luminosity is the observable most affected by 
the buildup of dense gas in cool cluster cores. The re- 
lations that are more sensitive to the ICM distribution 
on larger scales show less of an offset between CC and 
NCC clusters. In addition, these larger scale measure- 
ments (e.g., Micm(<?'50o)) show less scatter than rela- 
tions at smaller scales (e.g., AficM(<f'250o))- This is an 
interesting finding; it shows simply that core structural 
variations and temperature biases have a larger effect on 
relations measured in smaller regions around the core. 
This is an indication of the importance of core structure, 
which we discuss further below. 

2.3. K-hand Light-Temperature Relation 

The X-ray observables provide information about the 
ICM. We now turn to a very different cluster property, 
the K-hand galaxy light. The near- infrared (NIR) light 
traces stellar mass better than optical bands, and is well 
corre lated with such prope rties as cluster binding mass 
(e.g.. iLin et alJl20fl.'^ 1200^ . By comparing information 
from the X-ray and NIR observables, we hope to gain a 
better understanding of, for example, the true temper- 
ature bias introduced by cool cores. T he if-band data 
we use is taken from ILin et al.l l)2004f l: these data are 
available for 34 of the clusters in our sample. The mea- 
surements of near-IR light are effectively a simple sum of 
the light from individual galaxies, and thus represent a 
very different measure of galaxy cluster properties than 
the X-ray observables discussed above. 

At the bottom of Figure|21we plot the iiT-band luminos- 
ity projected within rspo- Here visual inspection suggests 
less evidence for a separation between the CC and NCC 
populations than in the X-ray relations. As with the X- 
ray observables, the scatter in the relation is dominated 
by the intrinsic scatter, as shown by the values for (Tint 
and (J law given in Table |^ 

3. COOL CORE AND NON-COOL CORE POPULATIONS 

It is clear that the presence of cool cores affects all 
of the ICM-related scaling relations. As previously dis- 
cussed, this is likely a result of various physical differ- 
ences between CC and NCC clusters, which include a 
simple temperature bias effect. Because we want to ex- 
amine merger-related scatter about the scaling relations, 
we must account for this separation of populations in a 
way that will allow us to compare CC and NCC clus- 
ters in a manner independent of cooling effects. One 
approach would be to treat these two populations indi- 
vidually, examining deviations from scaling relations in 
each. In this section we effectively take this approach 
by applying a correction to the CC population so that in 
the mean these clusters lie on the same scaling relation as 
the NCC population. Below we describe this approach as 
well as the amplitude of the offsets between the CC and 
NCC populations in each of the scaling relations. An al- 
ternative, more sophisticated approach will be presented 
in §12 

3.1. Aligning CC and NCC Populations 

To begin with, we align the CC and NCC population 
scaling relations by scaling the mean temperature of all 
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CC clusters by the same amount within a given relation. 
To find the appropriate scale factor for a relation, we in- 
crease Tx for the CC clusters by a range of factors (using 
the same factor for all CC clusters), re- measure the rel- 
evant observable at each temperature, and measure the 
reduced for the entire sample at each scale factor. 
Because the measurement uncertainties are so small, we 
use a combination of measurement and intrinsic scatter, 
which reduces the dependence of the scale factor on out- 
lier clusters. We measure versus temperature scale 
factor, then find a value for the intrinsic scatter CTint that 
makes the reduced = 1 at the minimum. We then 
find the temperature scale factor that corresponds to the 
new x^ minimum, and so on, until the process converges. 
That is, for an observable O, we find temperature scale 
factor A and intrinsic scatter dint such that 



TABLE 3 
Cool Core Temperature Scale 
Factors 



Scaling relation 



Tx scale factor 



-t'x(<»'2500)-7x 


1.47 




0.07 


ix(<»-50o)-Tx 


1.38 




0.06 


-'W'iCm(<»"2500)-Tx 


1.19 


± 


0.05 


-f^XCs(<''500)-Tx 


1.15 


± 


0.04 


^1.5xlO-13"'^X 


1.10 


± 


0.07 


MiCM(<r500)-7x 


1.05 


± 


0.04 


^3xlO-l*^'^X 


1.01 


± 


0.05 


-'^NIR(<f500)-7x 


0.77 


± 


0.08 



A^dof 



E 



\0,[\Ti)-0^,{\%)X 



= 1 



(4) 



at the x^ minimum for iVdof degrees of freedom. 
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Fig. 4. — X-ray luminosity projected within rsoo (ioy)^ and 
deviation of clusters from the best fit relation (bottom), after scal- 
ing the CC cluster temperatures to align the CC and NCC cluster 
populations. Open and filled markers correspond to CC and NCC 
clusters, respectively. The uncertainties include both the measure- 
ment uncertainty in the luminosity and an effective temperature 
contribution to the luminosity uncertainty as described in the text. 

Because aligning the CC and NCC cluster populations 
decreases the total scatter in each scaling relation, we 
must verify that the intrinsic scatter in the relations re- 
mains significantly greater than the measurement uncer- 
tainty. We again examine ix(<''50o)j after scaling the 
CC cluster temperatures by a factor of 1.38, determined 
by the method described. The resulting scaling relation 
is shown in Figure 01 (cf. the relation with non-scaled 
temperatures in Figure P). Although the total scatter 
in the relation is now significantly less than in the orig- 
inal relation, it still is clearly larger than the luminosity 
measurement uncertainties. We quantify this as before, 
measuring the raw and intrinsic scatter, and find that 
o'int= 0.33 and crraw= 0.34, again showing the dominance 
of the intrinsic scatter over measurement uncertainties 



(scatter measurements for CC temperature-scaled rela- 
tions are given in Table [JJ. 

A plot of x^ versus temperature scale factor for all 
of the X-ray and NIR scaling relations is shown in Fig- 
ure El the derived scale factors, plus their uncertainties, 
are listed in order from highest to lowest scale factor in 
Table |21 Clearly, a greater scale factor is required to 
align the CC and NCC populations for parameters that 
measure a smaller, more core-dominated region of the 
cluster, as discussed in § 12.21 Similar scale factors of 
~ 5% and ~ 1%, respectively, will align the two popula- 
tions for both MicM within r^oo and R^^iq-u. There is 
a somewhat larger difference between temperature scale 
factors for Micm within r25oo, ^i.sxio-i^i ^^'^ the core- 
subtracted luminosity within r5oo (~ 19%, ^ 10%, and 
~ 15%). A much greater scale factor is needed in ci- 
ther case for Lx 38% within r5oo, ~ 47% within 
''2500 )■ This demonstrates that the luminosity is affected 
by more than just the shift in temperature due to the 
presence of emission from a cool core. 



o 
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-0,4 -0,2 



Temperature shift factor 



Fig. 5. — Reduced versus CC temperature scale factor for 
each of the eight observable— temperature scaling relations. Pen- 
tagons: Projected luminosity within rsoo (filled) and r2500 (open). 
Squares: ICM mass within rsoo (filled) and r2500 (open). Tri- 
angles: Isophotal size for isophotes of 3 X 10 



arcmm ^ and 1.5 X 10 erg s ^ cm ^ arcmin (open) 



erg s ^ cm ^ 
2 ^'„„or^^ Stars: 



Projected core-subtracted luminosity within rsoo. Crosses: Pro- 
jected NIR light within rsoo- Vertical axis units are reduced 
value, but relations are offset vertically from one another for ease 
of viewing so vertical axis labels are not shown. 
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3.2. Examining the CC Temperature Bias 

Besides the cool core temperature bias, these scahng 
relations may be affected by cool core-related structural 
phenomena and merger-related phenomena that are in 
some way related to the presence or absence of cool cores. 
A cool, dense core results in a sharp central brightness 
peak that drives up the total luminosity of the cluster, 
adding to the separation between CC and NCC clusters 
on the Lx-Tx scaling relation. The increased central 
surface brightness also makes isophotal sizes larger and 
leads to higher measured gas masses in the central re- 
gions. However, idealized studies of mergers of spheri- 
cal, isolated clusters suggest that both ICM temperature 
and X-ray luminosity may be boosted during mergers, 
and that cluste rs can be perturbed from scaling relation 
by the merger ijRicker fc SaraziiJl200H) . If present, this 
behavior would also tend to separate the merging and 
non-merging (and hence NCC and CC) populations. 
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Fig. 6. 



1 5 10 

M,cM (<r5oo) [1013 MJ 

i?3xj^Q-i4 plotted versus A/icm(<'~50o)i along with 
the best-fit relation (top), and deviation in Rj from the best-fit 
relation (bottom). Open and filled markers correspond to CC and 
NCC clusters, respectively. Uncertainties are a combination of Rj 
measurement uncertainties and an effective uncertainty in Rj due 
to the mass uncertainty obtained using the scaling relation slope, 
as was done with temperature uncertainties for other relations. 

A good way of differentiating between these two effects 
is to examine a relation that has minimal temperature- 
dependence. To examine the impact of cool core effects 
on scaling relations at large radii, we plot R^xio-^"^ ^s^" 
sus Micm(<?'5oo) as shown in FigureEl Isophotal size is 
independent of temperature, and at a low isophote core 
effects should be of little importance. The ICM mass 
within a large radius will have only a slight dependence 
on core features, because only a small fraction of the 
cluster ICM mass lies within the core region. Though 
r^oo depends on temperature, any bias effect this intro- 
duces is weak. Hence, a scaling relation composed of 
these two observables provides a test of how much the 
cool core bias affects the structure outside the cluster 
core. Indeed, the relation shows no particular suggestion 
of separation between the CC and NCC populations, and 
the total scatter in the relation is much smaller than for 



any of the previously discussed observable-temperature 
relations. This provides evidence that the primary con- 
tributor to the CC/NCC separation in the observables at 
large radii is indeed a simple temperature bias, and not 
structural changes related to the development or disrup- 
tion of cool cores. This also gives us confidence that the 
offset between the CC and NCC populations is caused by 
the onset of a cooling instability within the cluster core 
rather than by shock-induced temperature and structural 
changes during mergers, which we would expect to be 
most apparent in observables that are sensitive to cluster 
structure outside the core (where relaxation timescales 
are the longest). 

Adopting this perspective, we can take the scale factors 
for relations involving observables that are less core sen- 
sitive to estimate the scale of the temperature bias. For 
example, the Micm(<''50o) and R^xio-^* relations have 
scale factors of 1.04±0.04 and 1.01±0.05, respectively, 
suggesting that temperature biases are at the few percent 
level, and that it is indeed structural differences in the 
core that are driving the larger scale factors seen in the 
more core sensitive observables. Interestingly, the near- 
IR relation shows less evidence for a cool core-related 
separation. In fact, following the same procedure as for 
the X-ray relations indicates that a negative scale factor 
of ^^23% is required to align the CC and NCC poplula- 
tions. (This may be partly driven by a few outliers, but 
having no reason to discard these data points, we do not 
exclude them.) This suggests that the galaxy population 
in CC clusters contains systematically less light than the 
galaxy population in NCC clusters, a result that deserves 
further attention. 

3.3. Intrinsic Scatter in CC and NCC Populations 

Having removed to first order, via temperature scaling 
of CC clusters, the separation between cluster popula- 
tions, we can begin to study the effects of merging on 
scatter about scaling relations; because mergers are ex- 
pected to disrupt cool cores, and because merging clus- 
ters are naively expected to have greater scatter about 
scaling relations, one would expect to observe greater 
scatter in NCC clusters. We therefore measure the in- 
trinsic scatter in the temperature-scaled CC population 
and the NCC population separately as a test of overall 
structural differences between them. That is, we mea- 
sure the scatter of each (CC and NCC) population about 
the same best-fit scaling relation. These values are given 
for all eight 0~Tx scaling relations in Tabled We use 
an F-test to quantify the significance of differences be- 
tween CC and NCC scatter for a given relation; the ta- 
ble lists the percent significance level at which equality 
of the variances is rejected. Remarkably, we do not ob- 
serve greater scatter in the NCC population; indeed, the 
CC population has significantly (i.e., > 68% significance) 
greater scatter than the NCC population in all but one 
scaling relation. 

The uniform CC temperature scaling method treats all 
CC clusters in exactly the same way. The greater scat- 
ter in the CC population, however, provides evidence of 
significant structural variation within that population — 
more variation, in fact, than in the NCC population. We 
are thus motivated to find a method to reduce CC/NCC 
separation that takes into account the variability of indi- 
vidual cluster structure. In the next section we present 
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TABLE 4 

Intrinsic Scatter CC and NCC Subsamples 



With CC Temperature Scaling 3 Paramter (O-Tx-Io) 

Scaling Relation CC aint — NCC (Ti„t — Uitf. CC (Ti„t — NCC (Ti„t — Uitf. 



Micm(<»-50o)-7x 


0.20 


0.09 


99.6 + 


0.17 


0.11 


93.0 + 


Micm(<?-2500)-Tx 


0.22 


0.10 


99.6 + 


0.17 


0.08 


99.7 + 


ix(<?-5oo)-'rx 


0.37 


0.28 


73.0 + 


0.28 


0.18 


90.7 + 




0.44 


0.31 


83.6 + 


0.40 


0.20 


99.3 + 


-^'XCs(<^500)-7x 


0.24 


0.24 


5.3 + 


0.22 


0.22 


3.6- 


^sxio-i'i-'^x 


0.15 


0.10 


91.9 + 


0.14 


0.11 


62.3 + 


-^1.5 X 10- 13 ""^X 


0.19 


0.10 


95.1 + 


0.16 


0.10 


84.0 + 


iNIR(<f'500)-Tx 


0.21 


0.09 


99.3 + 


0.23 


0.12 


97.3 + 



^Percent likelihood that scatter measurements for CC and NCC subsamples are different; see text. 
Plus sign indicates that CC scatter value is higher than NCC scatter value. 



such a method. 

4. PEAK SURFACE BRIGHTNESS AS A MEASURE OF 
COOL CORE STRENGTH 

In this section we discuss another method of reduc- 
ing cooling-related scatter in scaling relations: the use 
of peak surface brightness as an indication of cool core 
"strength" . We include surface brightness as a third 
parameter in scaling relations, show that this parame- 
ter's contribution to the relation is significant in all X- 
ray observable-temperature relations, and demonstrate 
its usefulness in reducing cool core-related scatter. 




10-13 10-12 lQ-11 

Iq [erg s"' cm"^ arcmin"^] 



Fig. 7. — Central cooling time from IPeres et all 1199811 plot- 
ted versus peak surface brightness Iq. Open and filled markers 
correspond to CC and NCC clusters, respectively. 

4.1. Brightness Measurements 

We use measurements of the peak surface brightness 
Iq from MME. These values were obtained by fitting (3 
models to azimuthally-averaged cluster surface bright- 
ness profiles. Clusters that appeared relaxed and dis- 
played residuals consistent with a central emission ex- 
cess were fit with a double (3 model, i.e., two models 
with the same f3 but different Iq and core radius were 
summed and fitted together; the cluster Iq is then the 
sum of the individual Iq from each model. In Figure |7| 



we plot the central cooling time versus Iq for the 45 clus- 
ters in our sample. There is a clear correlation between 
the two quantities, which is not surprising as the central 
cooling time is derived from the central surface bright- 
ness profile. Note that the PSF corrected central surface 
brightness varies by a factor of ^^500 for our flux limited 
cluster sample, suggesting that even in the low signal to 
noise regime it should be possible to differentiate the low 
and high surface brightness systems. 
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Fig. 8. — X-ray luminosity projected within rsoo plotted versus 
temperature (top), and deviation of clusters from the best fit Lx— 
Tx-/o relation calculated for three values of 7o {bottom). For this 
plot the cluster sample was divided into three subsamples based 
on Iq; the pentagons, squares, and triangles are the clusters in 
the subsamples with the highest, middle, and lowest values of Iq, 
respectively. Open and filled markers correspond to CC and NCC 
clusters, respectively. 

4.2. Observable-Temperature-Brightness Relations 

We now test whether Iq is a significant parameter by 
including it in the observable-temperature scaling rela- 
tions. That is, for each observable O we assume a func- 
tional form O oc T^I^ and solve for a, (3, and the nor- 
malization. We first examine the X-ray luminosity pro- 
jected within rsoo- A plot of ix(<r5oo) versus Tx is 
shown in Figure |S1 In this figure we have divided the 
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TABLE 5 

Best Fit Temperature and Brightness Scaling Parameters 



Observable Tx dependence Iq dependence Iq dependence 

(Tx bias removed)'' 



^ICm(<»'500) 


1.94 


± 


0.09 


0.06 


± 


0.02 


0.00 


± 


0.03 


Micm(<?'2500) 


1.92 


± 


0.08 


0.10 


± 


0.02 


0.04 


± 


0.03 


-I'X(<»'500) 


2.13 


± 


0.10 


0.26 


± 


0.03 


0.20 


± 


0.04 


-I'X(<''2500) 


2.33 


± 


0.18 


0.31 


± 


0.03 


0.24 


± 


0.04 


i'XCs(<''500) 


2.21 


± 


0.10 


0.10 


± 


0.03 


0.03 


± 


0.04 


^3x10-1^ 


1.03 


± 


0.07 


0.03 


± 


0.01 


0.00 


± 


0.01 


^1.5x10-13 


1.02 


± 


0.08 


0.07 


± 


0.02 


0.04 


± 


0.02 


LmR{<r60o) 


1.41 


± 


0.14 


0.00 


± 


0.03 


-0.04 


± 


0.03 



^That is, tlic factor 7 as defined in Eq.I^ 



cluster sample into three subsamples based on Jq, and 
we plot the measured Lx—Tx-Io relation for a value of 
Iq in the middle of each subsample; we also show the 
deviation in luminosity for each subsample. This plot 
shows qualitatively that the scatter about the Lx-Tx" 
Iq is much smaller than about the Lx"Tx relation (c.f. 
the Lx-Tx relation in Fig. Quantitatively, we find 
that the Lx{<r5oo)-Tx-Io relation has a best-fit power- 
law dependence on Iq with an exponent of 0.26 ± 0.03 
(uncertainty is obtained from bootstrap resampling and 
refitting); the Iq dependence is thus indeed significant. 
The raw and intrinsic scatter in ix(<f50o) about the re- 
lation are 0.26 and 0.24, respectively, much lower than 
the corresponding values for the original Lx-Tx relation 
and also lower than for the Lx-Tx relation with CC tem- 
perature shifting (see Table EJ. This dramatic reduction 
in scatter in the luminosity-related scaling relations sug- 
gests that it should be possible to use luminosity together 
with central surface brightness as a much more accurate 
cluster mass estimator than luminosity alone. 

O-Tx-Iq scaling relation exponents for the rest of the 
relations are given in Table |S1 All X-ray observables 
have a dependence on Iq that is significant at greater 
than two standard deviations; Lnir does not show any Iq 
dependence. We conclude that the peak surface bright- 
ness does provide us with useful information about clus- 
ter structure. Its usefulness for our present purpose is 
clear given the extent to which scatter about scaling re- 
lations is reduced by its introduction. Table[21gives mea- 
sured raw and intrinsic scatter values for each of the eight 
observable-Tx-/o relations, and separate measurements 
of CC and NCC scatter are given in Table |3| For X- 
ray observables we find reduced scatter in the O-Tx-Iq 
relations compared to the original O-Tx relations, and 
generally less than for the temperature-scaled relations. 
The exception is Lnir which, having no dependence on 
Iq, does not show reduced scatter compared to the orig- 
inal iNiR-Tx relation. 

There have been several recent studies of the evolution 
of cluster scaling relations at intermediate and high red- 
shift. Some studies have found positive evolution of clus- 
ter X-ray luminosity, consistent with that expecte d frorn 
self-similarity argum ents (e.g., IVikhlinin et al I i2002t 
IMaughan et alJl2005j) . but other observations are con- 
sistent with zero or somewhat negative evolution of Lx^ 
Tx and Micm"^^ scaling relation normalizations (e.g., 
I^rgani ct al. 2001; Holdcn et al. 2002; Ettori et all 
12004ft . Given the difference in normalizations for the CC 



and NCC populations, this may be partially explained 
by a simple change in the cool core fraction in the sam- 
ples being compared, which will shift the normalization 
of the entire cluster population. It would thus potentially 
be useful to compare O-Tx-Iq relations at low and high 
redshift, as the inclusion of a parameter measuring cool 
core strength may reduce this evolution effect. 

4.3. Temperature and Other Observable Biases 

We can use this three parameter scaling relation ap- 
proach to again estimate the scale of the temperature bi- 
ases. Scaling relations involving properties that have low 
dependence on core structure should have little to no de- 
pendence on Iq. That all X-ray observable-temperature 
relations show a significant Iq dependence can be taken 
as evidence for temperature biases in CC clusters; that is, 
the cool gas in CC cluster cores biases emission- weighted 
mean temperatures so that CC clusters appear to lie 
above scaling relations. We can attempt to quantify this 
bias in terms of Iq. We assume scaling relations of the 
form 

O = OqT^I(^ = Oq {\{Iq)TxT U (5) 
for each observable C If a scaling relation has no intrin- 
sic dependence on Iq (i.e., 7 = 0), then A(/o) oc 
We assume that the AficM(<''50o)~2x and R^ ^iq-h -Tx 
relations have no intrinsic Iq dependence (see § 13.211 , and 
so A(/o) oc /^o '''^='='^ oi)_ 'pj^g temperature bias thus varies 
by ^--^20% over the range of Iq in our sample; using the 
median values of Iq for the CC and NCC populations 
(4.5 X 10~^^ and 4.4 x 10~^^ erg s~^ cm~^ arcmin"^, 
respectively) gives an average CC temperature scale fac- 
tor of 1.07±0.02 (simply averaging the CC temperatures 
scale factors for Micm(<?'50o) and R-^y^iQ-u from 13.11 
gives 1.03 ±0.03). 

For the other scaling relations, 7 = /?— (0.03 ± 0.01)q. 
The values for 7, i.e., the dependence on Iq when the tem- 
perature bias is removed, are given in the fourth column 
of Table |31 Not surprisingly, the X-ray luminosities still 
have by far the strongest Iq dependence, which differs by 
a factor of four between the highest and lowest values of 
Iq in our sample. The Micm(<?'250o) dependence on Iq 
varies by ~30% over the sample. This suggests that the 
gas fraction varies significantly within r25oo depending 
on the strength of the cool core; variations of this scale 
will present challenge to those using cluster gas fractions 
of "relaxed" clusters to precisely study cosmology. 
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5. SUBSTRUCTURE AS A SOURCE OF SCATTER 

We now examine the relationship between substructure 
and the position of clusters on scaling relations. Having 
introduced two ways to remove the cool core tempera- 
ture bias, we now attempt to examine the merger related 
structural differences as deviations by clusters from scal- 
ing relations. In this section, we review the substructure 
measurements, and then discuss how deviation from scal- 
ing relations depends on substructure. We then compare 
results from the observational sample to an ensemble of 
hydrodynamical simulations. 

5.1. Substructure Measurements 

High-resolution instruments such as Chandra reveal 
hydrodynamic phenomena such as cold fronts that are 
clearly related to merging. However, it is generally not 
necessary to directly observe such features to find evi- 
dence of merger-related activity. Relatively crude, low- 
order moments of the X-ray surface brightness distribu- 
tion such as centroid variation and ellipticity have been 
shown to be effective at separating clusters with recent 
major me rgers from more relaxed systems (jjMohr ct ^ 
E993: Evr ard et al.l 1993), although these measurements 
are essentially unaffected by mergers along the line of 
sight. These measurements do not require the high res- 
olution of Chandra, and they were used to show that 
more than hal f of clusters display substructure in their 
ICM t Mohr et al. 1995). 

The centroid variation w is a measure of the "center 
shift" , or skewness of the photon distribution of a cluster. 
There are many ways to measure w; here, we measure 
within an isophote of 2 x 10~^^ erg s""'^ cm~^ arcmin"^, 
chosen as the lowest isophote that can be used for all 
of our cluster images. We measure the centroid of the 
portion of the cluster that has surface brightness above 
this isophote. We then examine the cluster at steadily 
brighter isophotes and measure the variance in the cen- 
troids measured for all these isophotes. The centroid 
variation w is the square root of this variance. We scale w 
to be in units of r^oo for each cluster; using the fractional 
variation rather than the raw variation simply accounts 
for the fact that the size of cluster virial radii can vary 
by as much as a factor of ^4 from low mass to high mass 
systems. Measuring the centroid variation using regions 
de fined by isophotes rather than using circular annuli as 
in iMohr et al.l l) 199,31) can provide information that cir- 
cular apertures may not; during mergers, clusters often 
do not have circular surface brightness distributions, and 
elongated structure may be missed or underemphasized 
by fixed circular apertures. 

We determine the axial ratio 77 from the flux-weighted 
second moments of the photon distribution, using an 
aperture centered on the brightness peak. That is, we 
measure moments 



(6) 



where the sum is carried out over all pixels within a cho- 
sen aperture, Xi are pixel coordinates {x or y) relative 
to the center of the aperture, and / is the measured in- 
tensity in the pixel. We measure 77 within an aperture 
of radius rsoo for each cluster, except for four clusters 
for which rsoo either is larger than our PSPC image or is 
close enough to the edge that background problems arise. 



In these cases we use apertures of radius r25oo- Using 
virial radii for the apertures provides a more physically 
meaningful scale for examining substructure than using 
a fixed metric radius. Diagonalizing the matrix obtained 
from equation @ gives the lengths of the major and 
minor axes, from which we then obtain the axial ratio. 
While high ellipticity is not a certain indicator of cluster 
substructure, and is sometimes observed even in appar- 
ently relaxed clusters (e.g.,|Hchueckcr et al. 2001), hydro- 
dynamical simulations show that during major mergers 
the ICM is typical ly highly flattened (.Evrard et al..il993 
IPearceet a,].lll994^. 
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Fig. 9. — Comparison of centroid variation w to the power ratio 

(vk) ipk) 

/Pq for the 45 clusters in the PSPC sample. Power ratio 
values are quoted in units of 10~^; centroid variation is plotted in 
units of r'soo- Open and filled markers are CC and NCC clusters, 
respectively. 

Another method to qua ntify cluster substr ucture 
makes use of "power ratios" (|Buote fc TsaillT995|) . This 
involves measuring moments of the surface brightness 
distribution E within some radius i?a 



I (-Rap) 



I](x')(i?')'"cosm(/)'d2a:' , (7) 



The "powers" Pm are then given by: 
Po = [aoln(-Rap)]' 



for TO = 0, and 



Prr 



2m2i?2m 



(8) 



(9) 



(10) 



for TO > 0, where i?ap is the radius of the circular aper- 
ture in which the moments are measured. The quantities 
of interest are the ratios Pm/Po; the division by Pq nor- 
malizes the flux within the radius of interest, allowing 
comparison of cluster observations with different fluxes 

and exposure times. The quantity p[^^'^ / P^^^\ which 
is calculated within an aperture centered on the cluster 
surface brightness peak, is similar to the centroid varia- 
tion w. Other ratios Pm>i/Po are measured within an 
aperture centered on the point where the centroid vari- 
ation is at a minimum; the quantity P2/P0 is related to 
the axial ratio rj. 
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Usin g a sample of PSPC observations, iBuote fc Tsail 
l)199fiD argue that certain relationships between power ra- 
tios may be viewed as evolutionary tracks. We seek here 
to find correlations between power ratios — mainly the 
simplest to interpret, p^P'^^ / p^P'''> and P2/P0 — and the 
deviations of clusters from scaling relations. We break 
from their approach of using fixed metric radii and use 
an aperture that scales with the cluster mass or temper- 
ature; specifically, we study the power ratios within the 
same characteristic radii (rsoo or r25oo) as we do with 
axial ratios. This provides a more physically meaning- 
ful scale for a cluster sample that spans more than an 
order of magnitude in mass. Figure |3 contains a plot 
of centroid variation versus P^^'^'' j p'^^^ and Figure ^] 
contains a plot of axial ratio versus PiIPq. The cen- 
troid variation w and the axial ratio r\ are correlated 
with the primary power ratio p'^^'^ jp''^^'^ and P2lPid^ 
respectively. Because neither pair of substructure mea- 
surements is perfectly correlated, we benefit from using 
all four measurements. 



axial ratio. It is clear that the Lx-Tx-Io relation is not 
masking any increase in scatter in high-substructure clus- 
ters. Similar results are seen for other scaling relations; 
we thus choose to focus on the O-Tx-Iq scaling relations 
for our study of merger-related scatter below. 
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Fig. 10. — Comparison of axial ratio r] to the power ratio P2/P0 
for the 45 clusters in the PSPC sample. Power ratio values are 
quoted in units of 10"'^. Note that axial ratio is plotted from high- 
est to lowest, i.e., (vertically) higher points have greater ellipticity. 
Open and filled markers are CC and NCC clusters, respectively. 



5.2. Substructure and Scaling Relations: CC 
Temperature Scaling versus O-Tx-Iq Relations 

To test for merger-related scatter in scaling relations, 
we measure the scatter about scaling relations by sub- 
samples of clusters, split according to the four substruc- 
ture measures discussed above. We wish to minimize cool 
core related scatter, and we have discussed two methods 
of doing so in this paper: uniform CC cluster tempera- 
ture scaling, and use of peak surface brightness as a third 
parameter in observable-temperature scaling relations. 

As shown above, the O-Tx-Iq relations generally have 
lower scatter than the CC temperature-scaled relations; 
however, one may wonder whether scatter information is 
being lost in the O-Tx-Iq relations. Qualitative com- 
parisons suggests that this is not the case. We give one 
example here: in Figure lTTI wc plot the difference between 
the data and best-fit ix(<^50o)~2x-/o relation and the 
best-fit CC temperature-scaled Lx{<r5oo)-Tx relation 
versus two substructure measures, centroid variation and 



Fig. 11. — Logarithmic (base-e) deviations of clusters from best- 
fit ix(<^50o)^^x^^O scaling relation (left) and best-fit I'x(<f 500)^ 
Tx relation (right) versus centroid variation w (top) and axial ratio 
rj [bottom). Open and filled markers are CC and NCC clusters, 
respectively. Centroid variations are given in units of rsoo . 
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Fig. 12. — Logarithmic (base-e) deviations of clusters from best 
fit scaling relations versus centroid variation w for each of the seven 
X-ray observable— temperature— brightness scaling relations and the 
NIR luminosity-temperature-brightness relation. Open and filled 
markers are CC and NCC clusters, respectively. Centroid varia- 
tions are given in units of rsoo- 

5.3. Substructure and Scaling Relations: Individual 
Cluster Relations 

We now look for merger-related structural variations 
in all scaling relations. Figure IT?! shows the natural log- 
arithmic deviation of each data point from the best fit 
scaling relation for our sample plotted versus centroid 
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TABLE 6 

Intrinsic Scatter in Scaling Relations, Split By Substructure 



Scaling Relation 








Split by w Split by r; Split by pf'^/p^^"'' 
L H Diff.(%)=' L H Diff.(%)^ L H Diff.(%)^ 


Split by P2/P0 
L H Diff.(%)=' 



Observational sample 



MicM(<r50o)-7x--fo 


0.18 


0.12 


95.5 + 


0.18 


0.13 


79.2 + 


0.15 


0.15 


2.5- 


0.15 


0.15 


8.4 + 




0.18 


0.11 


96.1 + 


0.15 


0.15 


4.6- 


0.14 


0.15 


17.9- 


0.14 


0.16 


48.9- 




0.32 


0.18 


98.4 + 


0.16 


0.31 


99.4- 


0.31 


0.17 


99.1 + 


0.30 


0.19 


94.9 + 


^-X ( <''2500 )-Tx-/o 


0.45 


0.19 


100.0 + 


0.32 


0.38 


57.8- 


0.38 


0.31 


61.7 + 


0.41 


0.28 


91.6 + 


i-XCS (<''50o)-'rx-/o 


0.26 


0.18 


86.2 + 


0.19 


0.24 


69.9- 


0.25 


0.19 


77.4 + 


0.24 


0.20 


57.7 + 


^3xlO-l*~'^X-^0 


0.16 


0.10 


93.1 + 


0.14 


0.12 


57.6 + 


0.12 


0.14 


51.9- 


0.14 


0.12 


42.4 + 


-^^LSxio-ia-Tx-Zo 


0.15 


0.14 


21.1 + 


0.09 


0.17 


98.3- 


0.13 


0.16 


56.6- 


0.13 


0.16 


65.4- 


i'NIR(<''500)-Tx-/0 


0.20 


0.21 


16.8- 


0.18 


0.22 


57.0- 


0.21 


0.20 


18.1 + 


0.19 


0.22 


43.9 + 


-'W'iCM (<''500 )--R3 X 10- 1-1 


0.06 


0.05 


37.1 + 


0.05 


0.06 


28.8- 


0.05 


0.06 


82.3- 


0.05 


0.06 


8.7- 


ixCs(<»'500)--R3xl0-l* 


0.05 


0.05 


10.5- 


0.05 


0.06 


28.2- 


0.05 


0.05 


37.6- 


0.05 


0.05 


15.1- 


-^^NIr(<^500)--R3x10-1* 


0.19 


0.14 


72.2 + 


0.20 


0.14 


82.0 + 


0.16 


0.17 


14.7- 


0.17 


0.16 


17.8 + 








Simulated cluster sample 














MicM(<r50o)-rx 


0.21 


0.21 


6.0 + 


0.21 


0.20 


10.9 + 


0.24 


0.17 


87.4 + 


0.22 


0.19 


39.5 + 


-f'XCs(<^500)-T'x 


0.26 


0.30 


56.1- 


0.25 


0.30 


59.3- 


0.23 


0.31 


78.1- 


0.27 


0.29 


29.4- 




0.09 


0.11 


62.6- 


0.10 


0.11 


54.8- 


0.11 


0.10 


12.4 + 


0.10 


0.11 


26.4- 



Note. — L and H arc low substructure (low w, P^^^^ / pI^^^\ and P2/-P0; and high 77) and high substructure, respectively. 



^Percent likelihood that scatter measurements for low and high substructure subsamples arc different; sec text. Plus sign indicates that 
low substructure sample scatter is higher; minus sign indicates that high substructure sample scatter is higher. 



variation, and Figure 1131 shows the deviation versus ax- 
ial ratio. The most obvious feature of these data is the 
semi-separation of CC and NCC clusters by substructure 
indicator; i.e., the CC clusters tend to have smaller cen- 
troid variations than the NCC clusters. This relationship 
between cool core status and axial ratio is not as striking, 
but is still present. 
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Fig. 13. — Logarithmic (base-e) deviations of clusters from 
best fit scaling relations versus axial ratio rj for each of the seven 
X-ray observable— temperature— brightness scaling relations and the 
NIR luminosity-temperature-brightness relation. Open and filled 
markers are CC and NCC clusters, respectively. Axial ratio in- 
creases to the left, i.e., ellipticity increases to the right. 



Figs . 1 1 2l and 1131 show no qualitative evidence for an in- 



crease in scatter in clusters with higher substructure, and 
indeed suggest greater scatter in clusters with less sub- 
structure. There may be a trend for the most irregular 
clusters to lie above the scaling relation, but the small 
number of clusters involved (2-3) makes this very un- 
certain. Cluster deviations from scaling relations versus 
p^(p'=)^p^(pfc)g^j-^j P^IPq (not plotted here) likewise show 
no suggestion of higher scatter in clusters with greater 
substructure. 

To quantitatively address the issue, we measure the 
intrinsic scatter in the scaling relations for subsamples 
grouped by centroid variation, axial ratio, P^^'^^ / P^^^\ 
and P2/Po- As before, we add a value for the intrin- 
sic scatter a-mt in quadrature to the uncertainty in the 
observable due to measurement and temperature uncer- 
tainty, and find the value of aint that results in a reduced 
value of unity for each scaling relation. We express 
(Tint in units of natural logarithm of the observable, i.e., 
the units of the vertical axes in Figs. 1121 and 1131 

We split the sample into two subsamples for each sub- 
structure measure; the split point for each is chosen to in- 
clude roughly half the clusters in the sample. Specifically, 
we split the sample at w = 0.02, rj = 0.875, p[P^'> / p^P''"' 
— 300, and P2/P0 — 20. Note that for the axial ratio 
77, a higher value corresponds to a more regular (spheri- 
cal) cluster, whereas for the other substructure measures 
a higher value is, roughly speaking, a messier cluster. 
Table contains the results; for ease of interpretation, 
a graphical representation of the same data is shown in 
Fig- El As in Table 01 we give the percent significance 
level at which equality of subsample variances is rejected. 

Broadly speaking, we find greater scatter in clusters 
with less substructure; with 11 relations and four sub- 
structure measures, we find significantly (i.e., same scat- 
ter rejected at > 68% level) greater scatter in low sub- 
structure clusters in 13 cases, and in high substructure 
clusters in only four. This is remarkable, although it is 
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not surprising in light of our earlier result (§13 that there 
is greater scaling relation scatter in the CC population 
than in the NCC population. One might suspect that 
a few outliers or a poor choice of splitting values could 
cause an apparent increase in scatter in lower substruc- 
ture clusters, but an examination of Figure IT?! does not 
support this. Indeed, while the scatter measurements 
suggest that clusters with higher ellipticity may have 
greater scatter. Figure [T^ suggests that this is a result of 
a few outliers, and does not constitute a general trend to 
higher scatter in more elliptical clusters. 
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Fig. 14. — Graphical representation of the observational sample 
data in Table |H] Open and closed bars are for low-substructure 
and high-substructure populations, respectively. 

One way to minimize the effects of the cool cores is to 
examine observables that are least affected by the cool 
cores. We examine deviation in Rj from the i?3x 10-14- 
-^icm(<''50o) scaling relation that was plotted in Fig- 
ure 1^ The emission- weighted mean temperature does 
not appear in this analysis, and these two observables 
are very insensitive to the core structure of the ICM. 
Deviations from this scaling relation are plotted ver- 
sus centroid variation in Figure 1151 Total scatter in 
this relation is smaller than for any of the observable- 
temperature relations, providing another indication that 
is it the core structure of clusters which varies most sig- 
nificantly within the population. The cluster deviation 
versus w provides some indication of higher scatter at 
lower substructure. We do see the same suggestion of a 
preferential boosting above the scaling relation for a few 
clusters at high w as for the observable-temperature rela- 
tions discussed above. We also examine scaling relations 
constructed from -R3xio-i<i-Lxcs(<'".50o) and i?3xio-i*" 
^nir(<?'50o); these, too, show little if any separation be- 



tween CC and NCC populations, and exhibit little evi- 
dence for different scatter for populations with different 
levels of substructure. 

In summary, the evidence clearly does not support our 
naive expectation that clusters with more substructure 
should exhibit higher scatter than their more relaxed 
counterparts. 
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Fig. 15. — As in Fig. 1121 but for deviation in Hj from the 
^3x10-i*^-'^iCm(<»'50o) scaling relation. 



5.4. Substructure and Multiple Scaling Relations 

Although we have shown that the clusters with the 
most substructure do not preferentially deviate from in- 
dividual scaling relations, it is possible that within the 
hyperspace defined by our broad range of observables 
these merging systems may have a tendency to lie some- 
what further from the general population. We examine 
this possibility by combining cluster deviations from all 
the scaling relations and probing for greater combined 
deviations in systems with the most substructure. To do 
this we assume that the cluster behavior about a scaling 
relation is a probabilistic indication of the consistency of 
that cluster with the typical structure of the population. 
Specifically, we assume that the probability of finding a 
cluster at a given deviation is described by a Gaussian 
in log space centered on the relation with a full-width at 
half-maximum equal to 2.354criogC), where criogo is the 
intrinsic scatter of the population about the scaling re- 
lation. That is, we calculate 
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for an observable O. We also calculate Ptotai, a measure 
of the significance of deviation from multiple scaling rela- 
tions together, by multiplying the individual probability 
amplitudes. If a cluster deviates slightly from three indi- 
vidual relations, for example, this should be reflected in 
.Ptotai- This approach assumes that the cluster behavior 
about each scaling relation is an independent indicator 
of its deviation from the whole population. 

We find no correlation between deviation from scal- 
ing relations and substructure for any individual O-Tx 
scaling relation, or for all observables measured together. 
The three highest substructure clusters noted in the pre- 
vious discussion do indeed have low total probability den- 
sities, as expected, but this is also true of several other 
clusters over the entire range of substructure in our clus- 
ter sample. We conclude that there is no readily dis- 
cernible relationship between substructure and deviation 
from scaling relations by individual clusters. 

5.5. Hydrodynamical Cluster Simulations 
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Although simulated clusters do not exhibit the full 
complexity of real clusters, carrying out our analysis on 
a controlled sample of well-understood systems is an im- 
portant component of our work. We use a simulated clus- 
ter ensemble consisting of 45 smoothed particle hydrody- 
namics simulations evolved in a ACDM cosmology; de- 
tails of the tech niques used and of this pa rticular sample 
can be found in lBialek et"an l)2001'.'200^. These simula- 
tions have resolution suffici ent to exhibit cluster merger 
features such as cold fronts ( Bial ek et al.i2 0021 . They do 
not include any ICM cooling mechanism, and hence cool 
cores will be absent; this provides a good opportunity to 
examine the results that would be expected in our obser- 
vational sample if we were able to completely remove the 
cool core effects. Lack of resolution and incomplete phys- 
ical modeling make the central, core regions of simulated 
clusters unreliable, and so we choose to examine scaling 
relations that are less sensitive to these core regions. 



-^^iCM(<f'50o) using the actual simulation data (i.e., not 
calculating Micm from mock observations). We mea- 
sure the isophotal size in the ROS AT hand correspond- 
ing to an isophote of 1 x 10~^ counts s'-^ arcmin"^; 
this instrumental isophote leads to isophotal sizes that 
approximately match the normalization of the observed 
^3xio-i*~^x relation. This isophote generally lies well 
outside the core of the simulated clusters. 

Scaling relations for these three observables are shown 
in Figure 1161 As we did for the observed scaling rela- 
tions, in Figure [T7I we plot the difference between each 
cluster and the best-fit relations versus the w and 77 sub- 
structure indicators. There is some suggestion here of a 
trend toward greater scatter at higher substructure. To 
quantify this, we calculate the intrinsic scatter as was 
done for the observations, both for the entire sample and 
for subsamples of roughly equal size; values of CTint are 
shown in Table 
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Fig. 16. — Simulated scaling relations for AficM within 
^"200 1 core-subtracted X-ray luminosity projected within rsoo, and 



isophotal size for an isophote of 1 X 10 
versus Tx. 
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We examine the projected ROSAT-hand (0.5-2.0 keV) 
core-subtracted luminosity within rsoo; that is, without 
the luminosity projected within 0.20r5oo, as in our ob- 
served cluster sample. We derive the scaling relation for 
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Fig. 17. — Deviation in base-e logarithm of simulated clusters 
and from best-fit relations plotted versus centroid variation w and 
axial ratio r] for each of the three simulated X-ray scaling relations. 
Centroid variation is given in units of rsoo. 

Looking at all four substructure measures, there is a 
tendency toward lower scatter in clusters with less sub- 
structure. However, in general the quantitative differ- 
ences are rather small; significantly (> 68% significance) 
different subsample scatter is found in only two relations, 
and greater scatter in low substructure clusters in only 
one. We also calculate the probability amplitude for indi- 
vidual clusters by combining information from all scaling 
relations, as we did for the observations. This approach 
provides no suggestion of greater or smaller deviation for 
clusters with less substructure. Thus, while the simula- 
tions do not show evidence for greater scatter in clusters 
with more substructure, they also do not show the ten- 
dency toward higher scatter in less messy clusters that 
we see in the observed cluster sample. 

The scatter we measure in the simulated sample can 
be compared directly to the intrinsic scatter in the ob- 
served samples. Because of the lack of radiative cooling 
in the simulations, it is likely best to compare the sim- 
ulation scatter to the cool core "corrected" observations 
(i.e., the scatter in the O-Tx-Iq relations). Agreement 
in the scatter of simulated and observed scaling relations 
would serve as one more indication (along with the slope 
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and normalizations: iBialek et al.ll200l(l that the simula- 
tions are an accurate representation of real clusters. In- 
terestingly, the scatter is greater for the simulations in 
mass (0.20:0.14) and luminosity (0.27:0.21), and smaller 
in isophotal size (0.10:0.13). The comparison is not en- 
tirely appropriate in the case of the ICM mass, because in 
the simulations we use the actual three dimensional ICM 
mass measured within rsoo, and for the observations we 
calculate this through a deprojection. The bottom line is 
that although these simulations with early and uniform 
preheating do a good job of reproducing the observed 
slopes and amplitudes of scaling relations, there is still 
clearly work to be done to match the scatter in observed 
scaling relations. The range of missing physics in the sim- 
ulations includes radiative cooling, conduction, magnetic 
fields, feedback from AGN, galaxy formation, etc., and 
all of these must be considered together with the effects of 
finite spatial resolution. It is quite interesting to consider 
that some ingredients currently missing from simulations 
could actually reduce the scaling relation scatter — that 
is, reduce the cluster to cluster structural variations at a 
given mass. 

5.6. Summary of Substructure Results 

We have shown that clusters with greater substructure 
do not preferentially lie significantly farther from scaling 
relations than clusters with less substructure. In fact, 
in observed clusters, there is a tendency toward greater 
scatter about scaling relations in apparently more re- 
laxed clusters. These findings contradict the naive ex- 
pectation that cluster structure is greatly disturbed by 
merger events and so should increase scaling relation 
scatter. 

Comparison of the observational and simulated cluster 
samples suggests that the greater scatter in more appar- 
ently relaxed clusters must be the result of processes not 
present in the simulations. Cool core-related phenomena 
are clearly the prime candidate for the higher scatter seen 
in observed clusters with less substructure, especially as 
our results show unquestionably higher scatter in CC 
clusters when compared to NCC clusters. AGN activ- 
ity may also contribute, as AGNs occur frequently in the 
cluster population and can produce radio cavities with 
assoc iated energies of at least ~ 10^° erg l|Birzan et al.l 
l2004|) . 

However, even the simulations show only weak evi- 
dence of higher scatter in clusters with more substruc- 
ture. We must conclude that either there are mecha- 
nisms which introduce a range of structural variations 
into apparently relaxed clusters, or that mergers sim- 
ply do not perturb cluster structure to the extent ex- 
pected and suggested by simulations of isolated clusters 
ijRicker fc Sarazinll200l[|Randall et alJl200^ . 

It may be suggested that the use of archival temper- 
atures from different sources may introduce scatter or 
otherwise hide merger or cooling effects. Redoing the 
ix(<'"50o) portion of our analysis with 43 of 45 temper- 
atures taken from a single published source iWhitell200(H) 
demonstrates that this is not the case; while quantitative 
scatter measurements can be sensitive to outliers, we find 
no evidence that the qualitative trends we report are af- 
fected by moderate shifts in cluster temperatures. 

6. CONCLUSIONS 



We examine the relationship between substructure and 
cool cores in galaxy clusters and the scatter about X-ray 
and near-IR scaling relations using emission-weighted, 
non-spatially resolved, non-cooling corrected tempera- 
tures. We separate clusters into CC and NCC subsam- 
ples according to their central cooling times, and we 
quantify substructure using the centroid variation, ax- 
ial ratio, and power ratios p'^^^'^ j p^^^^ and P2I Pq- With 
these tools and a wide range of crude cluster observables 
and uncertainties, we examine a large number of galaxy 
cluster scaling relations in an approximately X-ray flux 
limited sample of 45 clusters. Our primary findings are: 

1. There is an offset between CC and NCC popula- 
tions on all observable-temperature scaling rela- 
tions that we examine. This separation is partly 
due to an emission weighted mean temperature bias 
at around the 7% level. The offset must also be due 
to differences in the core structure of the two sub- 
samples, because those scaling relations that are 
most sensitive to the cluster core show the largest 
offsets. This offset is not driven by recent merg- 
ers, because scaling relations involving observables 
that are sensitive primarily to the outer structure 
of clusters show larger scatter in the CC popula- 
tion. 

2. We show that the central X-ray surface brightness 
can be used to characterize the "strength" of cool 
cores, and that introducing it as a third parameter 
in observable-temperature scaling relations greatly 
reduces the scatter about those relations. Thus, the 
central surface brightness provides a tool for study- 
ing the evolution of cluster scaling relations in a 
manner less sensitive to any change in the fraction 
of cool core clusters with redshift. In addition, the 
small scaling relation scatter when using the sur- 
face brightness means that crude cluster observ- 
ables like the X-ray luminosity and temperature 
can provide more accurate virial mass estimates 
than are obtained without the third parameter. 

3. Parameterized in terms of central surface bright- 
ness /q, we find that the emission weighted mean 
temperature bias correction factor is proportional 
to /^''■03±0-0i)^ Given the factor of ~500 variation 
in /o within our sample, this implies a maximal 
differential correction across our sample of ~ 20%. 
Using the median /q for the CC and NCC popula- 
tions, we calculate an average CC temperature bias 
factor of 1.07 ± 0.02 for this sample. 

4. We find that although CC clusters tend to exhibit 
less morphological substructure, they exhibit at 
least as much scatter as (and often more than) the 
NCC clusters. Thus, structural variations among 
CC clusters are at least as large as the structural 
variations among NCC clusters. This result has 
important implications for analyses which rely on 
the presence of a cool core to indicate that a cluster 
is relaxed. 

5. Clusters with greater morphological substructure 
do not exhibit more scatter about scaling relations 
than clusters with lower substructure. In fact, we 
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observe a trend toward higher scatter in clusters 
with less substructure. This may be partially due 
to the structural variations within cool cores, which 
are typically found in clusters that exhibit less sub- 
structure; however, even after using central surface 
brightness to reduce the systematic cool core ef- 
fects, we find that the clusters with less substruc- 
ture exhibit as much or more scatter as those with 
more substructure. 

6. The differences between low and high morphologi- 
cal substructure clusters are modest in our study of 
hydrodynamical cluster simulations without cool- 
ing; there is only a weak indication that those clus- 
ters with higher substructure exhibit higher scat- 
ter. 

7. As in the purely X-ray scaling relations, there is 
no relationship between deviation from the Lnir- 
Tx relation and cluster substructure. However, in 
the NIR relation there is a negative temperature 
scale factor required to align the CC and NCC 
populations, and a lack of any dependence on Iq 
in a constructed iNiR~2x-/o relation. Because we 
know that emission-weighted mean temperatures 
are biased by the cool core gas in CC clusters, this 
suggests a difference between galaxy populations in 
CC and NCC clusters. A detailed study of differ- 
ences in NIR properties of the galaxy population in 
CC and NCC clusters would clearly be interesting. 

Together, these results from studies of real and simu- 
lated clusters indicate that cool core related phenomena 
(such as radiative cooling, AGN activity, and entropy 
injection at an earlier epoch), and not cluster merging, 
are the primary sources of scatter in scaling relations. 
Perhaps it should come as no surprise that X-ray ob- 
servables that arise from emission, which is sensitive to 
the square of ICM density, are most perturbed by the 
detailed structure of the cluster core. However, the lack 
of a strong substructure related enhancement of scat- 
ter in scaling relations without sensitivity to the cluster 
core and in hydrodynamical simulations of clusters, is 
surprising. It suggests that perhaps all clusters retain 
departures from equilibrium at a significant enough level 
that even recent mergers do not perturb their structure 
enough to make them appear unusual. A young popu- 
lation with relaxation timescales that are comparable to 
the time since the last major merger would presumably 
exhibit this kind of behavior. 

Reconciling these observations and our conclusions 
with high res olution hydrodynamical mergers of idealized 
clusters fe.g.. lRicker fc Sarazinll20r)]l) requires a rarity of 
such large-scale mergers or perhaps additional physics 
within the ICM that suppresses the boosts. Note that 
correlated excursions in luminosity and temperature dur- 
ing a merger will not suffice as an explanation, because 
we have shown using 8 different observable-temperature 
scaling relations that there is no strong relationship be- 
tween substructure and scaling relation scatter. We have 
demonstrated here that lower resolution hydrodynami- 
cal simulations simply do not produce large deviations 
from the general cluster population even when there is 
evidence for significant substructure, either in individual 



clusters or in the high-substructure cluster population 
as a whole. Whatever the explanation, it is clear from 
reasonably large samples of real clusters that there are 
no outliers on the scale of those predicted by the high 
resolution, idealized cluster merger simulations. 

Our result may be consistent with the scenario where 
NCC clusters evolve to beco me CC clusters in the ab- 
sence of major mergers fe.g.. lOta et al.ll2005|) . However, 
the larger scatter we observe for CC clusters in all scaling 
relations raises important questions. In particular, even 
after using the central surface brightness to correct for 
CC effects, we still observe higher scatter about the CC 
relations (see Table0}. Moreover, we see the larger scat- 
ter in CC clusters even in scaling relations that involve 
observables that are not sensitive to the core structure 
(i.e., faint isophotal size, ICM mass within rsoo, and NIR 
light). Thus, if CC clusters evolve from NCC clusters 
because of an absence of mergers, then the observations 
require some other source of scatter or variation in clus- 
ter structure to be present throughout the cluster virial 
region. 

Alternatively, cool cores may arise through a scenario 
that is driven by something othe r than th e recen t merger 
history of the cluster. McCart hv et al.l l)2004|) suggest 
that variations in entropy injection into the intraclus- 
ter medium could determine whether or not a cool core 
forms. This varied entropy injection would also con- 
tribute to structural variations or scatter in scaling re- 
lations. Because our cluster sample indicates that it is 
the CC clusters which exhibit the highest scatter around 
scaling relations (even those relations with little core sen- 
sitivity), it seems likely that this or some similar, non- 
merger driven scenario is responsible for the presence or 
absence of cool cores in clusters. Within this scenario 
the tendency for cool core clusters to exhibit less mor- 
phological substructure would be primarly due to the 
effects of the often dominant X-ray bright core, which 
would tend to bias axial ratios high and centroid varia- 
tions low. That is, in cool core clusters a morphological 
substructure indicator is in large part reflecting the char- 
acteristics of the bright, symmetric core. 

It will be quite interesting to return to this question of 
the dual nature of galaxy clusters — youthful as indicated 
by the high frequency of morphological substructure, yet 
strikingly regular as indicated by scaling relations — with 
new tools and larger samples extending over a wider 
range of redshift. Of particular interest will be the addi- 
tional leverage afforded by the new generation of high 
signal to noise Sunyaev-Zel'dovich effect observations, 
which should be dramatically less core-sensitive than X- 
ray observations. With the tens of thousands of clusters 
expected in dedicated surveys, it should be possible to 
quantify with high significance any subtle, merger related 
trends that may be present. 
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